Diffusion-limited deposition with dipolar interactions: fractal dimension and 

multifractal structure 



in 
o 
o 

(N 
o 

CD 

Q 

o 



M. Tasinkevych,^^^ J. M. Tavares,^'"* and F. de los Santos''' 

^ Max- Planck- Institut fiir Metallforschung, Heisenbergstr. 3, D-70569 Stuttgart, Germany 
^Institut fiir Theoretische und Angewandte Physik, 
Umversitdt Stuttgart - Pfaffenwalding 57, 70569 Stuttgart, Germany 
'^Centra de Fisica Tedrica e Computacional da Universidade de Lisboa 
Avenida Professor Gama Pinto 2, P-1649-003 Lisbon, Portugal 
^Instituto Superior de Engenharia de Lisboa 
Rua Conselheiro Emidio Navarro, 1, P-1949-014 Lisboa, Portugal 
^ Departamento de Electromagnetismo y Fisica de la Materia, 
Umversidad de Granada, Fuentenueva, 18071 Granada, Spain 
(Dated: February 2, 2008) 

Computer simulations are used to generate two-dimensional diffusion-limited deposits of dipoles. 
The structure of these deposits is analyzed by measuring some global quantities: the density of the 
deposit and the lateral correlation function at a given height, the mean height of the upper surface 
for a given number of deposited particles and the interfacial width at a given height. Evidences 
are given that the fractal dimension of the deposits remains constant as the deposition proceeds, 
independently of the dipolar strength. These same deposits are used to obtain the growth probability 
measure through Monte Carlo techniques. It is found that the distribution of growth probabilities 
obeys multifractal scaling, i.e. it can be analyzed in terms of its /(q) multifractal spectrum. For 
low dipolar strengths, the f{ct) spectrum is similar to that of diffusion-limited aggregation. Our 
results suggest that for increasing dipolar strength both the minimal local growth exponent Qmm 
and the information dimension D\ decrease, while the fractal dimension remains the same. 
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I. INTRODUCTION 

The formation of clusters and deposits by irreversible 
aggregation of particles is an example of a nonequilib- 
rium growth process. Although several mechanisms can 
be involved in these processes, the most simple mod- 
els only take into account the effects of thermal diffu- 
sion. Among these models, the DLA (diffusion limited 
aggregation)^ and DLD (diffusion limited deposition)^'^ 
have been widely used, fn DLA, particles are released one 
by one at a distance R from a seed particle, and perform 
a random- walk in a d-dimensional space. Eventually, the 
random walker either reaches the aggregate and attaches 
to it, or moves sufficiently far away from the aggregate to 
be removed. DLD is a version of DLA for the growth of 
deposits on fibers and surfaces. ^•'^ In DLD, particles also 
diffuse randomly through a rf-dimensional space, but at- 
tach either to a d — 1 dimensional substrate or to the 
deposit. Despite being based on a very simple algorithm, 
these models (and modifications thereof) generate very 
complex patterns and have served as a guideline for the 
understanding of wide range of phenomena such as elec- 
trochemical deposition, viscous fingering, dendritic solid- 
ification, dielectric breakdown, etc.^ 

At present, there is no theoretical framework to de- 
scribe the scaling behavior of DLA or DLD structures. 
Initially, it was assumed that the structure of the ag- 
gregates is homogeneous, statistically self-similar fractal 
and could be characterized by a fractal dimensionality 
D. However, this simplified assumption does not fully 
capture the complexity of DLA and DLD structures and 
a better characterization is obtained from the studies of 



the growth probability distribution. One can then apply 
a multifractal scaling analysis to determine the multifrac- 
tal spectrum of the growth measure, or equivalently an 
infinite hierarchy of fractal dimensionsj^ The remark- 
able scaling behavior of the growth probabilities^ imme- 
diately raises the question of their behavior when mech- 
anisms other than the thermal diffusion arc included. 

It is well known that short-range isotropic interac- 
tions do not change the cluster's structure at large length 
scales. At small length scales, however, a short-range at- 
traction (repulsion) promotes the formation of less dense 
(more dense) aggregates but with no change in the frac- 
tal dimensions D.^'"^'^ A different picture emerges when 
long-range, anisotropic interactions are considered. Re- 
sults for DLA of dipolar particles for d = 2 indicate that 
D decreases from D Ri 1.7 (the value for pure DLA^) to 
Z3 « I.l when the dipolar interaction is increased.^" This 
is in accordance with the results for cluster-cluster aggre- 
gation of dipolar particlesii as well as with experimental 
results for the aggregation of magnetic micro-spheresii^ 

In our previous work^— we have grown two-dimensional 
DLD clusters consisting of up to 10^ dipolar particles, 
and found evidences of a more complex behavior for the 
deposits' structure. Fig. ^ shows two typical deposits 
obtained for weak and strong dipolar interactions. The 
deposits consist of many tree-like clusters competing to 
grow. As the deposition process continues fewer and 
fewer trees keep on growing. Eventually only a single 
tree survives. In both cases the height and width of an 
individual tree of size s (number of particles in the clus- 
ter) can be described by the power-laws H ~ s"^" and 
W ^ s'^^ , respectively^ Likewise, the average number 
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FIG. 1: Typical DLD deposits which are obtained for a) weak 
and b) strong dipolar interactions. The deposits contain a) 
12000 and b) 2000 particles. 



Us of trees of a given size scales with s asris ^ s^'^ }^ Two 
scaling regimes have been observed: (i) for s less than the 
crossover size s*, the shape and the fractal dimension of 
the trees Dt are temperature dependent in such a way 
that T and Dt decrease with increasing of the interaction 
strength. We call this scaling regime the dipolar regime. 
(ii) For s > s* , and for large enough systems, pure DLD 
values of the exponents are observed independently of 
the interaction strength. This implies a crossover to the 
diffusion-driven DLA scaling regime, where the effects of 
the dipolar interactions are dominated by thermal effects. 
The dipolar regime corresponds to deposits exhibiting an 
orientational order of dipoles and the onset of the DLA- 
regime coincides with the disappearance of the orienta- 
tional order. More strikingly, it has been found that the 
value of D (the fractal dimension of the entire deposit) 
barely varies with the dipolc interaction strength. In the 
present paper, we will show further evidence of this by 
analyzing in detail the scaling behavior of the density and 
a lateral correlation function. 

Despite having the same fractal dimension, deposits 
obtained for different strengths of the dipolar interac- 
tion (equivalently, different effective temperatures) ex- 
hibit quite different structures. This situation resembles 
that of DLA and percolating clusters in d = 3, namely 
they have the same fractal dimensions, although their ge- 
ometrical structures are quite different piSiii In this paper 
we characterize the deposits of dipoles further by testing 
a multifractal scaling of the growth probabilities {pi}, 
where pi is the probability that the perimeter site i is 
the next to be added to a cluster. 

This paper is organized as follows: in section^ we in- 
troduce the model of dipolar DLD and describe details 
of the simulations. In section IIII Al the scaling expo- 
nents obtained for the density, lateral correlation func- 
tion, the mean height of the upper surface and the in- 
terfacial width are compared with those of DLD. In sec- 
tion ^^B] we present the multifractal spectra calculated 
for strong and weak dipolar interactions, and for several 
stages of growth. Finally, in section IIVI we summarize 



our findings. 



II. MODEL AND SIMULATIONS DETAILS 

The model and the simulation technique are the same 
as in our previous worksiii^; we briefly outline them 
here, and describe simulation details which allowed us to 
grow relatively large clusters. The simulations are per- 
formed on a two-dimensional square-lattice with lattice 
spacing a and width La. The adsorbing substrate corre- 
sponds to the bottom row {y = 0) and periodic boundary 
conditions are applied in the direction parallel to the sub- 
strate (x-direction). Particles carry a three-dimensional 
dipole moment of strength /i and interact with other par- 
ticles through the dipolar pair-potential (pD'- 

2) = — 3- [3(/ii • fi2)(/i2 • ^"12) - Ml ■ A2] : (1) 
'"12 

where ri2 = \fi - > a, fi2 = (n - ■r2)/ri2, n is a 
two-dimensional vector and p,i, (12 are three-dimensional 
unit vectors in the direction of the dipole moments of 
particles 1 and 2 respectively. 

Initially, a particle with random dipolar moment is in- 
troduced at the lattice site (a;™, Hmax + AL), where Xin 
is a random integer in the interval Hmax is the 

maximum height of the deposit and ^ is a constant. The 
particle then diffuses by a series of jumps to nearest- 
neighbor lattice sites, while interacting with the particles 
that are already part of the deposit. Eventually, the par- 
ticle cither contacts the deposit (i.e., becomes a nearest- 
neighbor of another particle that belongs to the deposit), 
or attaches to the substrate (i.e., reaches the bottom of 
the simulation box), or moves away from the substrate. 
In the latter case, the particle is removed when it reaches 
a distance from the substrate larger than Hmax + "^AL, 
and a new one is introduced. Once the particle attaches 
to the substrate or to the deposit, its dipole relaxes along 
the direction of the local field created by all other parti- 
cles in the deposit. In all simulations reported here we 
take A = 1; larger values of A have also been tested 
and found to give the same results, but with drastically 
increased computational times. 

The effect of the dipolar interaction on the random 
walk is incorporated through a Metropolis algorithm. 
The interaction energy of an incoming particle with the 
deposit is given by E{M,f,(i) = J2iii4'D{i, M + 1), 
where r denotes the particle's position, p, is the orien- 
tation of its dipole moment and M is the number of par- 
ticles in the deposit. To simplify the notation, the sum 
over the periodic replicas of the system is omitted in the 
last expression. To move the particle, we select a neigh- 
boring site r and a new dipole orientation p,' randomly. 
This displacement is accepted with probability 

_ . \ E{M,f',p')-E{M,r,p)]\ 
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where T* = ksTa^/fi^ is an effective temperature in- 
versely proportional to dipolar energy scale. In the limit 
T* 0, only displacements which lower the energy E 
are accepted. On the other extreme, for T* — > oo, all dis- 
placements are accepted and our model reduces to DLD. 

One can estimate the number of calculations of 
t, necessary to grow a deposit of mass M on a strip of 
width L in the following way. Since each incoming parti- 
cle starts its movement at a distance of order L from the 
deposit, it takes roughly steps to reach the deposit. If 
there are n particles in the deposit, has to be calcu- 
lated nL^ times. Summing this factor over M particles 
gives t ~ L^M^. Since we have performed simulations for 
L ~ 10^ and M ~ 10^, these values roughly give t ^ 10^^, 
which is about 3 to 4 orders of magnitude larger than the 
corresponding values of the largest systems simulated in 
an equilibrium run of two-dimensional dipolar fluids»i^ 
Therefore, to make simulations feasible we rewrite the 
expression for the incoming particle-deposit interaction 
energy in a different form. To this aim, the dipolar pair 
potential Eq. is rewritten as = jiDij-fij, 

where Dij is the dipolar field created by the particle i at 
the site occupied by the particle j 

Dij ^-A- [3(A» • r^j)r^j - Aj] , (3) 

Using Eq. Q one can write the interaction energy be- 
tween the incoming particle (which is at site fj and has a 
dipole orientation fij ) and a deposit formed by M dipoles 

as 

£;(M,f,,A,) = M[>(M,j)-A,, (4) 

where U{M, j) = ^ij (^g^™j the sum over peri- 

odic replicas is omitted) is the dipolar field created by 
the deposit at site fj. During simulations we store the 
dipolar field U{j) at site rj together with the size of the 
deposit M{j) for which the field U is calculated. The 
particles in the deposit are ordered according to their ar- 
rival "time" . Then, the interaction energy is rewritten in 
the form E{M, rj,fij) = n{U{j) + 'EZmU)+i ^^o) ' Aj, 
which allows to reduce the number of calculations of the 
dipolar pair-potential from M to M — M{j). At the early 
stages of growth most of the lattice sites have not yet 
been visited by an incoming particle, hence M{j) — 0. 
However, as the deposit grows the number of the lattice 
sites that have been visited more than once increases. 
Moreover, due to the shadowing effect, the sites close to 
the tips of the trees become visited more frequently, and 
thus the overall computation time of the dipolar energies 
is decreased. We have performed tests to compare the 
computation times for calculating the energy between 
the two methods. For a system of width L ^ 10"^ and 
M ~ 10'^ particles, the storage of the dipolar fields rep- 
resents gains in computational time of about 3 orders of 
magnitude. 

Finally, the long-range of the dipolar interaction is 



treated by the Ewald sum method adapted to the slab 
geometry of the systemii^ 



III. RESULTS 

Simulations were carried out using four temperatures, 
T* = 10-1, 10-2, 10-3, lO-'', and four system sizes, L = 
200,400,800, 1600 with 20 000, 30 000, 50 000, and 100 
000 particles per deposit, respectively. Each choice of 
these parameters corresponds to one of the two regimes 
of growth. For instance, the DLA regime is never ob- 
served for the lowest temperature T* = 10-^, while for 
T* = IQ-i it is easily reached even for the smallest sys- 
tem size L = 200. As we mentioned above, our pre- 
vious findings suggest that the fractal dimension of the 
deposits D is the same as for DLD even in the dipolar 
regime^ To validate this picture, we next examine the 
scaling behavior of the following quantities: the density 
of the deposit and the lateral correlation function at a 
given height; the mean height of the upper surface for a 
given number of particles; the interfacial width at a given 
height. All of them are global quantities that characterize 
the entire deposit, and which do not show any qualitative 
or quantitative differences between the dipolar and the 
DLA regimes. 



A. The fractal dimension 

A typical plot of the mean density p{h) at a distance 
h from the substrate has three regimes separated by two 
crossover heights, hi and hs (see Fig.EJ. At early "times" 
the deposit builds up until it reaches a height hi. Then, 
there appears a scaling regime during which the density 
decreases as a power of height, p{h) ~ h~°' , the exponent 
a being related to the fractal dimension D of the deposit 
by D = 2 — The density stops decreasing and sat- 
urates when the lateral correlation length ^|| reaches the 
size of the system L. Given that ^|| grows with height as 
^11 ~ /i^, then hs ^ L"^ , with 7 = C This behavior can 
be described by the scaling form 

pih,L)^L-^f{h/L-'), (5) 

where f{x) is a scaling function with the properties: 
f{x) ^ for small x and f{x) —^ const for large x. It 
is clear that (3 = aj. A linear regression between hi and 
hg for the largest system-size (L = 1600) gives roughly 
the same scaling exponent 0.25 < a < 0.29, showing no 
systematic variation with temperature. This implies a 
value of the fractal dimension 1.71 < D < 1.75, which 
is the same as for DLD. As we mentioned above, this is 
in contrast with a previous study, Ref. 0, where a con- 
tinuous variation of -D as a function of T* was reported. 
Although several source of discrepancy are possiblepSi we 
believe that the estimates of Ref. suffer from insuffi- 
cient statistics and rather small cluster masses, that pre- 
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FIG. 2: Mean density p of deposits as a function of the dis- 
tance h from the substrate. The densities are calculated for 
a system-size L — 1600 and several values of the effective 
temperature T*. Beyond the crossover height hs, the density 
saturates to a constant value pa. The horizontal and vertical 
dashed-lines indicate the ps and hs obtained for T* = 10~*. 




FIG. 3: Lateral density- density correlation function calcu- 
lated according to Eq. © for several distances h from the 
substrate. The system-size is L = 1600 and the effective tem- 
perature T* — 0.1. The position of the minimum Xmin of the 
correlation function as a function of h is shown in the inset. 
a;min scales with h with a scaling exponent ^ ~ 0.82. 



vent the true asymptotic regime (M ~* oo) from being 
reached. 

According to Eq. ||SJ), for ft, » L the mean density 
attend a constant value ps that scales with the system- 
size as Ps ^ L^^. For DLD the exponent /3 « 0.33i^ 
This value agrees with those obtained for T* = 10"'' and 
T* = 10~^ (the only temperatures for which density sat- 
uration is neatly observed), /3 = 0.33(1) and fi = 0.34(1), 
respectively. From the scaling relation ^ — P/a one has 
7 = 1.3(1). An independent check of this value was car- 
ried out by calculating the height hs at which the den- 
sity saturates and which is related to the system-size by 
hs ^ L'^ . There is substantial uncertainty in the determi- 
nation of hs and, again, we are limited to the two lowest 
temperatures, T* = lO""* and 10~^. For both of them 



7 = 1.30(5), while for DLD a value of 7 « 1.20 has been 
measured. 

There is a third route to the exponent 7 through 
the two-point lateral density-density correlation function 
C{x, h) for horizontal cuts at a height h. This correlation 
function is defined as 



C(a;,ft) = -^^a(a;' 



x, h)a{x' , ft). 



(6) 



where (j{x, h) is an occupation number which equals 
unity if the lattice site {x, ft) belongs to the deposit and 
zero otherwise. Fig. |31 depicts a plot of C{x, h) as a func- 
tion of X for several heights ft. The main feature is a 



pronounced minimum at Xr, 



ft^ which can be in- 



terpreted as the mean distance between the trees at the 
height ft, Ref. IT3. and that, for DLD, scales with ft with 
an exponent C, ~ 0.80 — 0.85. For very large heights, the 
exponent Q increases slowly with increasing height and 
may approach unityiii Owing to poor statistics, our sim- 
ulation data enables us to evaluate C only for T* — 10~^ 
and 10~^, thereby extending our estimations of 7 to the 
entire range of available temperatures. Our results are 
consistent with C, ~ 0.85, equivalently 7 k, 1.2, and no 
significant deviation from DLD behavior is found. 

Some other aspects of the morphology of the deposits 
have also been investigated. In our earlier work^ the 
mean height ft of the upper surface h{x) (h{x) is defined 
as the maximum height of the occupied sites which are in 
the column x) was shown to scale with the number of par- 
ticles M (in the large-M limit and for every temperature) 
with an effective exponent (f> — l/{D—d+l) « 1.33—1.44, 
which corresponds to a fractal dimension D 1.69 — 1.75. 
In the present paper we double the number of particles 
of the deposits to 10^, which allows an improvement of 
these bounds to D !v 1.71 — 1.73. As we wrote in Ref.lTM 
had the deposits been allowed to grow only to interme- 
diate stages (e.g., lO"* particles) an apparent variation of 
D with T* would have obtained. 

Owing to scale invariance, the exponent (f) should coin- 
cide with the growth exponent which controls the power- 
law divergence of the width of the upper surface W = 
([ft(a;) - ft]2)i/2 _ M't'. Linear fits yield « 1.3 - 1.4, in 
agreement both with our previous estimation and with 
the value we have measured for DLD, </) = 1.34(2). The 
divergence of the width is cut off at the the system-size 
length scale, and the width reaches a saturation value 
Wsat ~ , where the exponent i5 is a conventional mea- 
sure of interface roughness. Unfortunately, for every tem- 
perature saturation is barely reached only for L — 200, 
whilst for the rest of sizes much bigger deposits would be 
necessary. 



B. Multifractal spectrum in the dipolar and the 
DLA scaling regimes 

We discuss hereafter the dependence of the multifractal 
spectrum of the growth-measure on the system size, the 
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total number of particles, and the effective temperature. ability measure 



1. General remarks 

For DLA and related models, the growth probabil- 
ity varies sharply from site to site, changing from large 
values at the outer tips to very low values for the less 
accessible points deep within the cluster. It has been 
found that long-range repulsions enhance the screening 
effects, which leads to the formation of less branched 
clusters. Contrarily, long-range attractions reduce the 
growth probability at the cluster's tips and lead to the 
formation of more compact objects. In Fig. 21a we show 
the iso-lines of the visit probability distribution in the 
vicinity of a deposit grown at T* = 10"''. To estimate the 
visit probabilities we use probe particles that follow bi- 
ased random-walk trajectories, and register the frequency 
of visiting the unoccupied lattice sites. To illustrate the 
effects of the dipolar interaction, in Fig.0|b we show the 
visit probabilities around the same deposit estimated un- 
der purely diffusive conditions (T* = cx)), i.e. the dipolar 
interaction between the probe particles and the parti- 
cles in the deposit has been turned off. The interaction 
strongly enhances screening behavior of the cluster. More 
perimeter sites become screened, i.e. they are character- 
ized by smaller values of the sticking probability. On the 
other hand, the dipolar interaction enhances the concen- 
tration of the growth-measure at the outer sites of the 
cluster's branches. 

Studies of the growth probability distribution for 
DLASSiSSiSi and related modelsS^ indicate that the distri- 
bution of the growth-measure can be described in terms 
of a multifractal scaling model. Suppose that the bound- 
ary of a cluster of size L is covered by S{e) boxes of size e. 
Then, a growth measure p{i) is introduced in the i'th box 
as the probability for a random- walker to land in box i. 
As mentioned above, p(i) varies sharply over the cluster 
surface in such a way that the local growth probability 
density diverges at the tips and goes to zero inside the 
"fjords" (both behaviors are cut off at the particle length 
scale). To each box one can associate a singularity expo- 
nent a{i) via 



(7) 



In what follows we keep e fixed at the lower cut-off length 
scale, which is of one lattice unit e = 1, and L/e is varied 
by examining systems with different sizes L. The be- 
havior of the number of boxes N{a)da with the scaling 
exponent a taking on a value in the interval [a, a -I- da] 
defines the local scaling density exponent /(a) 



N{a) - L^("\ 



(8) 



One then introduces the scaling function T{q) which char- 
acterizes the scaling behavior of the moments of the prob- 



i 

Using Eqs. iQ and ©, Z{q) can be written in the form 
Z{q) - / daN{a)p'^ = I dai^("'-«". (10) 



Evaluating this integral by saddle-point approximation 
leads to 



Z{q) - i/("(9))"9°('3)^ 



(11) 



where the functions a{q) and f(a{q)) are defined implic- 
itly by the condition 



dfja) 
da 



(12) 



By comparing Eq. (|ll|l with Eq. the following relation 
is obtained 



-Tiq) = f{a{q)) ~ qa{q)- 



(13) 



The moment scaling function T{q) is related to the fam- 
ily of dimensionalities Dq introduced by Hentschel and 
Procaccia^^ through T{q) — {q ~ 1)-Dq. The limit Dq = 
D = limq^o+ Dq is the fractal dimension of the clus- 



Dn is the information dimen- 



ter, and Di = lim^_^i+ j^q 
sion. For DLA, T{q) is a non-linear function of q, i.e. 
an infinite hierarchy of exponents is required to charac- 
terize the moments of the probability measure. Some 
points of the f{a) spectrum can be related directly to 
the Dq: the maximum of / is given by /(a;(g = 0)) = -D, 
f{a{q =1)) = a = Di, a,„„i = D^c, and a^ax = £'-oo- 

Usually, T is calculated as a function of q and then the 
multifractal spectrum f{a) is obtained after performing 
a Legendre transform. The multifractal spectra that we 
present here, however, are obtained through interpola- 
tion of numerical histograms, using equations |(7J and 

. Because of the unknown normalization constants in 
these scaling relations, this method provides local growth 
exponent a and local scaling density /(a) up to additive 
corrections of order ~ 1/lnL, what causes the / spectra 
to fail to satisfy a number of important properties in- 
cluding data collapse, tangency of the f = a line and the 
f{a) curve at a single point (corresponding to g = 1 and 
giving l?i), or mislocation of other representative points. 
Here, we shall show that taking into account these correc- 
tion terms improves considerably the measured spectra 
(at least for low a). 

With this purpose, first note that for a finite L all 
quantities p, N,a, f,T must carry a label L. We define 
the size dependent ul and /l as following: pL = L^"^ , 
Nl{oil) = L^^("^\ To bring in the corrections to scaling, 
the integral H10|) is evaluated retaining the next-order 
terms to give 



Zdq) 



£-^-^■(9) K LfLiaL)-qaLiq)- ""^JP -G{L) 



(14) 
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where G(L) = ln[lnL/(27r)]/(21nL) and fl stands for 
the second derivative of with respect to a^. Notice 
that 



daL 



(15) 



holds. Assuming now that, to order 0(1/ Ini), tl k, t 
and using Eqs. 114|l . (|15|l the foUowing expressions can be 
readily derived 



flit 

J L 



2(/i')MnL' 



(16) 




3.0E-09 2.0E-07 1.9E-06 



/(a)=/L(aL) + 



fi fill 
Jl-Il 



H-fD 

2\nL 



G{L). (17) 



Taking as the independent variable, instead of 9, then 
these two equations give the parametric representation of 
the asymptotic multifractal spectrum /(a) as a function 
of measured quantities /l , q;l . 



2. Calculation of fi^{aL) trough numerical histograms 



FIG. 4: Iso-lines of the visit probability distribution in the 
vicinity of a deposit of dipoles. The probabilities are esti- 
mated numerically by using probe particles and counting the 
frequency of visiting of each unoccupied lattice site. The de- 
posit contains 2000 dipoles and has been grown at T* = 10~*. 
In a) the strength of the dipolar interaction of the probe parti- 
cles with the dipoles in the deposit corresponds to an effective 
temperature T* — 10~*. In b) the interaction is turned off 
(T* = 00) and the probe particles follow random-walk tra- 
jectories. Dark and light regions correspond to high and low 
visit probabilities, respectively. 



The growth probabilities PLii) are estimated numeri- 
cally using probe particles that sample the growth mea- 
sure. The probe particles follow biased random- walk tra- 
jectories in the dipolar field created by the cluster. Then 
p{i) is calculated as PL{i) = NL{i)/NT, where iVi(i) is 
the number of trajectories which have terminated on the 
perimeter site i and Nt is the total number of trajectories 
(probe particles). Having estimated the growth proba- 
bilities PL{i), we then calculate the number of perimeter 
sites iVi(ai)AaL with the value of Phii) in the interval 
where is defined as 



lnpL(0 
InL ' 



(18) 



More specifically, we calculate the histogram A^l(q;l) 
from which the multifractal spectrum is obtained as 



/i(aL) 



\nNL{aL) 
InL 



(19) 



In the calculations we choose for the width of the bin 
Aaz, ~ 0.15. 

High effective-temperatures results are shown in Fig.|Sl 
/l(q:l) is obtained at T* = 0.1 for deposits of size 
L = 400, 800, 1600 and number of particles M in the 
range lO'^ < M < 10^. The number of probe parti- 
cles used varies from 10^ for {L,M) = (1600,10^) to 
1.6 X 10^ for (L,M) = (400,3 x 10^). The collapse of 
the curves shows evidence that the distribution of the 
measure is multifractal, with /l(q:l) gradually converg- 
ing towards a well-defined spectrum as the deposit mass 
increases. An illustration that the uppermost curves are 
asymptotic is shown in Fig. \^ which displays data only 



for deposits of size L — 400 and containing M ~ 2000, 
10000, 20000, 30000 and 50000 particles (notice the small 
change in when AI changes from 30000 to 50000). 
The minimal value of the local growth exponent aL is 
associated with the region where the measure is most 
concentrated. It is found that aLmin ~ 0.56, a value 
that compares favorably with measurements for DLA for 
d — 2, a„iin « 0.59(4)1^ amin is related to the fractal 
dimension of the deposit D through D > 1 + amin^ 
This inequality is in agreement with the values of D es- 
timated by the scaling of the densities. One can identify 
a region of small ckl's for which the multifractal spec- 
trum does not change with the deposit mass M. In par- 
ticular, the value of aLmin ~ 0.56 is the same for the 
biggest and the smallest deposits, (L = 1600, M = 10^) 
and (L = 400, M = lO'^), what suggests that the maxi- 
mum growth probability, within the considered range of 
M, does not depend on M. Beyond the boundary value 
~ 1.1 there is clear mass dependence. Finally, no- 
tice that fL{ctL) flattens for larger values of a^. The 
maximum possible value of /l, fmax, should be equal 
the fractal dimension D of the deposit. For the DLA 
the maximum has been reported to be at the position 
ao ~ 4.5, Ref. With the Monte-Carlo technique we 
can not determine such small probabilities. However, one 
can anticipate that fmax is not far from the expected 
value D. 

Low effective temperatures results are displayed in 
Fig. El fL{aL) is calculated at T* = lO'^ for both the 
initial and late stages of growth (small and large M), 
the number of probe particles varying from 7.9 x 10^ 
for (L,M) = (1600,5 x 10^) to 5 x 10^ for (L, M) = 
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FIG. 5: fLioLL) spectrum of the growth probabiUty mea- 
sure for DLD of dipolar particles calculated at an effec- 
tive temperature T* = 0.1 as defined in section This 
value of T* corresponds to the DLA scaling regime for all 
values of M presented in the figure. ctL ~ — ln(pL)/lnL 
with PL being the growth measure, (a) Results obtained 
for system-sizes L = 400, 800, 800, and deposit masses in 
the range 10^ < M < lO'^. (b) demonstrates the conver- 
gence of fhiotL) spectrum towards an asymptotic regime for 
M > 2 X 10**. We find a minimal value of the local growth 
exponent ULmin ~ 0.56. The solid straight lines are given by 
/l = O-L- 



(400,3 X IC*). The initial stages of growth correspond 
to the dipolar regime, with the deposits containing 
M = 1000,2000,4000 particles for L = 400,800,1600, 
respectively. At the late stages of growth {AI — 
30000,50000,50000 for L = 400,800,1600 respectively), 
the DLA regime has been attained, as described in sec- 
tion IIIA. The inset of Fig. demonstrates the conver- 
gence of fhiciL) towards the asymptotic behavior which 
is reached at M > 20000. Again, like in the case 
T* — 0.1, a region of small ol's can be identified where 
the multifractal spectra does not depend on M . In fact, 
in both regimes we obtain a value of aLmin ~ 0.36, which 
is considerably lower than the corresponding value ob- 
tained for T* ~ 0.1. The turning point separating M- 
dependent from Af -independent behavior is located at 
aL ~ 1.35. Beyond this point good data collapse is still 
observed for the dipolar regime, whereas the collapse of 



FIG. 6: /L(a_t) spectrum of the growth probability measure 
for DLD of dipolar particles calculated at T* = 0.001, with 
QfL ~ — ln(pL)/lnL. The results are obtained for system sizes 
L = 400, 800, 1600, and deposit masses in the range 10^ < 
M < 5 X 10". The solid curve schematically represents the 
spectrum for T* = 0.1. The inset demonstrates the 
convergence of the fL(Q-L) spectrum towards an asymptotic 
regime for A/ > 2 X 10*. The values of deposit masses M < 
4000 correspond to the dipolar scaling regime; the individual 
trees of the deposits have the fractal dimension Dt ~ 1.33. 
For deposit masses 2 x 10* < Af < 5 x 10* the system is 
in the DLA scaling regime with Dt ~ 1.56. We observe the 
minimal value of the local growth exponent aLmin ~ 0.36 in 
both dipolar and DLA scaling regimes. 

the curves is not impressive after the crossover to the 
DLA regime has taken place. We shall see below how 
to amend this shortcoming. The changes in the shape 
of the /i-spectrum as Af grows can be interpreted as an 
increase in the fractal dimension of the screened parts of 
the deposit, but for the studied range of Ul's, /l still 
stays well below the T* =0.1 curve. According to the 
picture developed in section Till AI the maximum values 
of /l should equal D « 1.71. The figure suggests a fur- 
ther increase in if more probe particles are used to es- 
timate the growth probabilities. Although no conclusive 
statements can be made regarding this point, we shall see 
in the next subsection that it is expected that the large 
parts of the spectra shift upwards once corrections to 
scaling have been included. 

Lastly, the elongated appearance of the trees in the 
dipolar regime resembles that observed in the initial 
stages of deposits grown at high temperatures in the DLA 
regime. Our data show, however, that the multifractal 
structure of the dipolar regime is clearly different from 
those of the initial stages of growth of deposits at higher 
temperatures, even if the fractal dimensions are similar. 



3. Corrected spectra 

When the multifractal measure possesses a continuous 
spectrum, the straight line ?; = a is tangent to the f{a) 
curve at f{a) = aM> This general behavior is not seen 
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FIG. 7: Corrected spectra f{a) calculated according to the 
Eqs. (113,1113 for T* = 0.1, upper panel, and T* = 0.001, 
lower panel. The numerical curves fhiaL) are fitted to the 
function ao + ai exp(a2a^^) + 04 ln(ai:,). The symbols repre- 
sents the fitted functions, and the lines - corrected spectra 
/(q). For T* — 0.1 the numerical data only in the range 
ctLmin < OLL < 2.1, aud for T* = 0.001 in aLmin <aL< 1.5, 
are fitted. 



in the hitherto shown spectra because they have been 
obtained through numerical histograms. Hence, the scal- 
ing relations expressed by Eqs. 0, provide the local 
growth exponent a and the local scaling density f{a) up 
to additive constants ~ 1/lnL. This is the reason why 
the curves shown in the Figs. [Sj |H| instead of being tan- 
gent to the line y = a, i ntersect it. The multifractal 
spectra presented in Refs. Il7l25l29l were also calculated 
as numerical histograms, and display the same behavior. 

From a more physical point of view, /(a) = a is a 
turning point that separates regions of the spectrum cor- 
responding to high (small a) and low (large a) growth 
probabilities. Moreover, the subset satisfying /(a) = a 
has a fractal dimension equal to the information di- 
mension Di of the multifractal formalism introduced by 
Hentschel and Procaccia.^^ The particles in this subset 
carry almost all the growth probability, in such a way 
that Di is the fractal dimension of the active zone, i.e. 
the unscreened surface. The identification of Di from 
Figs. 13 is troublesome since there is no tangency at 
Jl = ctL- The use of the corrected expressions (|16ll and 



(|17|l is therefore called for, but their application immedi- 
ately poses the problem of computing up to third order 
numerical derivatives from noisy, scarce simulation data. 
We therefore resorted to fitting to the simulation data 
several functional forms. In Fig. [7| results for the cor- 
rected f{a) are shown for the particular choice of the fit- 
ting function aQ + ai exp(a2Q;2^) + lii(aL). The numer- 
ical data are fitted only for the values of in the range 

aLrmn < CtL < 2.1 for T* = 0.1, and aLmin < CkL < 1.5 

for T* — 0.001. For the larger values of ul there is a big 
uncertainty in the numerical data, thus those parts of 
the fLipii,) curves are excluded from this particular fit- 
ting procedure. We stress that the Fig. [3 is shown only 
as a qualitative illustration of the effects of the correc- 
tions terms in equations (|15|l and H17|) . In particular, the 
tangency property is restored in every case, while data 
collapse is preserved. Estimating Di as the value of / 
where the /-curve has unit slope we obtain Di « 0.98 
for T = 0.1, which is close to the exact result Di — 1 
for DLAi2£ With decreasing T* we observe the decrease 
of Di, namely we find Di w 0.84 and Di « 0.70 for 
T* = 0.01 (the corresponding figure is not shown) and 
T* = 0.001 respectively. Other functional forms yield 
very similar results. 

The corrected spectra have other interesting features 
as compared to the uncorrected ones. At this point, how- 
ever, the discussion has to be kept at a qualitative level 
since at the extreme values of one faces two different 
numerical problems: i) on approaching amin the deriva- 
tives of any order of "true" spectra / diverge; ii) on the 
other hand, large ul values correspond to data with poor 
statistics. Some general observations can nevertheless be 
made. There is strong indication that, upon incorpora- 
tion of the corrections terms: i) the part of the /^-curve 
corresponding to large aL shifts upwards for every tem- 
perature; ii) data collapse emerges for T* — 0.001. In 
particular, for the values of a w 3.0, /(a) takes simi- 
lar values for both values of T*, a result consistent with 
the interpretation that the fractal dimension does not de- 
pend on the effective temperature. Note also that a re- 
markable improvement in data-collapse can be observed 
for the case T* = 0.001 and large M. Finally, there is 
a slight tendency for amin to move to the right. For 
T* = 0.1 the corrected value gets even closer to its DLA 
value whereas for T* — 0.001 it keeps at amin ~ 0.4. We 
stress that on the whole this same behavior carries over 
to other fitting choices. We also calculate the corrections 
to the measured /L(aL) spectrum in the dipolar scaling 
regime (T* = 0.001) and small values of M). As ex- 
pected we get the same amin and Di as for large M, but 
no concluding evidence that D is going to be the same. 

According to the Turkevich-Scher conjecture D = 
1 -|- Umin,'^^ such bchavior of amin would be in contradic- 
tion with the fact that D does not change with effective 
temperature (as concluded in section Fill A|l . The above 
mentioned relation between D and amin is obtained us- 
ing implicitly the assumption that the most extremal 
sites of the clusters, or tips, are the most active ones. 
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i.e. atip = Oirain- We have checked the location of the 
sites with maximum growth probability in the clusters of 
simulations, and found that they often do not coincide 
with the clusters tips. Thus, aup > am.in and the more 
general relation D > 1 + amin , Ref • 1271 applies in this 
case. 



IV. CONCLUSIONS 

The diffusion-limited deposition of magnetic particles 
shows a crossover between two regimes, with a crossover 
size that sensitively depends on the temperature. At the 
early stages of growth both the fractal dimension of the 
trees and their size-distribution, as given by the ex- 
ponent r, are temperature dependent and have signif- 
icantly smaller values as compared with those of pure 
DLD (dipolar regime). As the size of the trees exceeds 
the crossover value s* the diffusion-driven DLA scaling 
regime emerges. In this regime Dt and r have the same 
values as in DLDjH Here, we have provided evidences 
that the fractal dimension of the entire deposit remains 
constant as the deposition proceeds. This has been done 
by analyzing the density profile of the deposit, the lateral 
correlation function, and the mean height and the width 
of the upper surface. It ensues that Dt and r conspire so 
as to give a fixed value of D. 

Multifractal analysis shows that, for each value of the 
interaction parameter and at each stage of growth, the 
normalized distribution of growth probabilities can be 
scaled onto a single curve using the same scaling form 
as in DLD thereby providing the /(a)-spectrum. The 
features of the f{a) spectra allows us to distinguish the 
structure of the deposits at high and low temperature (see 
Fig. P). We have found that amin decreases significantly 
with decreasing temperature, revealing that the concen- 
tration of the growth probability becomes more and more 
marked when dipolar interactions are increased. Like- 
wise, the fractal dimension of the active zone Di de- 
creases with decreasing temperature, meaning that for 



low temperatures less sites are involved in the growth of 
the deposit. As a consequence, and since Di and amin 
do not depend on the stage of growth, the presence of 
dipolar interactions reveals in the structure of the de- 
posits by increasing further the probability of growth at 
"hotter" sites and by originating less dense deposits. On 
the other hand, the f{a) spectra at late stages of growth 
suggest that the fractal dimension of the deposits does 
not depend on temperature (or dipolar interaction), cor- 
roborating our previous results. However, our findings 
also indicate that the /(a)-spectra in the high- and low- 
temperature DLA scaling regimes are different. 

The multifractal spectrum in the dipolar regime (low 
temperature at the early stages of growth) is clearly dif- 
ferent from that in the initial stages of growth at high 
temperature, thus providing evidence that these two sit- 
uations have to be held distinct. The f{a) spectrum 
obtained for the dipolar regime was, however, not accu- 
rate enough to confirm that the fractal dimension of the 
deposits in this regime is approximately equal to that of 
later stages of growth. 

In a more general perspective, this work shows that 
the information dimension Di and the scaling exponent 
a^nin are much easier to determine through a numerical 
measurement of the f{a) spectrum then the fractal di- 
mension D, provided that the corrections of Eqs. H16() 
and (|17|) are taken into account. In fact, the low a part 
of a spectrum can be obtained with good statistics using 
small systems (see Figs. El and EJ. Therefore, our work 
suggests that the effect of interactions in DLA (or DLD) 
deposits could be more easily studied through Di and 

a^iin • 
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